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Abstract 

We report on a first next-to-next-to-leading order calculation of the decay constants of the D (D*) 
and B (B*) mesons using a covariant formulation of chiral perturbation theory. It is shown that, using 
the state-of-the-art lattice QCD results on /d s //d as input, one can predict quantitatively the ratios of 
Id*/ Id*, IbJIb, and Jb*/ Ib* taking into account heavy-quark spin-flavor symmetry breaking effects on 
the relevant low-energy constants. The predicted relations between these ratios, /dj/ Id* < /d s //d and 
/b s //b > fDs/ Id, and their light-quark mass dependence should be testable in future lattice QCD simula- 
tions, providing a stringent test of our understanding of heavy quark spin-flavor symmetry, chiral symmetry 
and their breaking patterns. 

PACS numbers: 12.39.Fe, 13.20.Fc, 13.20.He, 12.38.Gc 
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The decay constants of the ground-state D (D*) and B (B*) mesons have been subjects of 
intensive study over the past two decades. Assuming exact isospin symmetry, there are eight 
independent heavy-light (HL) decay constants: f D (f D *), f Ds (f D *), f B (f B *), f Ba (f B *). In the 
static limit of infinitely heavy charm (bottom) quarks, the vector and pseudoscalar D (B) meson 
decay constants become degenerate, and in the chiral limit of massless up, down and strange 
quarks, the strange and non-strange D (B) meson decay constants become degenerate. In the real 
world, both limits are only approximately realized and, as a result, the degeneracy disappears. 

The gluonic sector of Quantum ChromoDynamics (QCD) is flavor blind, so the non-degeneracy 
between the HL decay constants must be entirely due to finite values of the quark masses in 
their hierarchy. A systematic way of studying the effects of finite quark masses is the heavy- 
meson chiral perturbation theory (HM ChPT) \ The HL decay constants have been calculated 
up to next-to-leading order (NLO) in the chiral expansion, and to leading-order (LO) |4j, |5Q and 
NLO in 1/ nriH expansion, where m B is the generic mass of the HL systems. In a recent work, 
a covariant formulation of ChPT has been employed to study the pseudoscalar decay constants up 
to NNLO for the first time and faster convergence compared to HM ChPT was observed [|8|]. 

Lattice QCD (LQCD) provides an ab initio method for calculating the HL decay constants. 
There exist many rtj = 2 + 1 computations of the pseudoscalar decay constants, f Dg and fn [|9- 
131. and f Bs and f B [Q,[l4 .[l^l, motivated by the important role they play in determinations of the 
CKM matrix elements and in tests of the standard model (see, e.g., Ref. Illql ). On the other hand, 
for the vector meson decay constants, most existing simulations are quenched Ill7l - ll9ll . except for 
Ref. [20] where nj = 2. Simulations with nj = 2 + 1 are underway [21]. 



In this letter, we report on a first next-to-next-to-leading order (NNLO) covariant ChPT study 
of the HL pseudoscalar and vector meson decay constants. We will show that heavy-quark spin- 
flavor symmetry breaking effects only lead to small deviations of the ratios f B Jf B , /d*//d*, 
and f B *J f B *, from f D j ' f D . Utilizing the latest HPQCD data on f Ds and f D [10], and taking 
into account heavy-quark spin-flavor symmetry breaking corrections to the relevant low-energy 
constants (LECs), we are able to make some highly nontrivial predictions on the other three ratios. 
The predicted light-quark mass dependencies of the HL decay constants are also of great value for 
future lattice simulations. 

The decay constants of the D and D* mesons with quark content qc, with q = u,d, s, are 
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defined as 



(0\q^ l5 c(0)\P q (p)) 
(0\q^c(0)\P*(p,e)) 



-if P y, 



(1) 

(2) 



where P q denotes a pseudoscalar meson and P* a vector meson. In this convention, fp has mass 
dimension one and Fp* has mass dimension two [12211 . For the sake of comparison with other 
approaches, we introduce fp* = F P */m P *, which has mass dimension one. Our formalism can be 
trivially extended to the B meson decay constants and therefore in the following we concentrate 
on the D mesons. 

To construct the relevant Lagrangians in a compact manner, one introduces the following fields 1 
and currents as in Ref. PO: 

i$ + mp 



H 
J 



2m i 



-(Yp; + tPi b ), 



where P = (D°, D+, D+), P* 
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r)*+\ J I juc jdc Tsc\T 
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with the weak current 



JT = 97/i(l — 7 5 ) c ' m P i s me characteristic mass of the P triplet introduced to conserve heavy 
quark spin-flavor symmetry in the mp — > oo limit: mp, at NLO and mc at NNLO (see Table 

Mu^d^u + ud^v)) and 



r„ with r„ 



1). The covariant derivative is defined as 
u 2 = U = exp[||] with $ the pseudoscalar octet matrix of Nambu-Goldstone (NG) boson fields, 
F their decay constant in the chiral limit. The weak couplings have the following form yfl: 



r(2) 

r(3) 



aTr[J b H a ]ul b , 

£ [i^Ti[J b H4 ab ] + ^Tr[J b d„H a ]u: b 
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mp 
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{6 D Tr[J 6 # a ](x + « t ) a6 
+6 A Tr[J 6 iy a ]^ 6 (x H 



where a is a normalization constant of mass dimension two, u h 



(5) 
(6) 

(7) 

AttFo is the 



scale of spontaneous chiral symmetry breaking, and Xh 



M 



&\&^{jr^,m\,'lm 2 K — m 2 ). Here and in the following Tr denotes trace for the Dirac matrices. 



1 It should be noted that the heavy-light states in the relativistic formalism have mass dimension of 1 instead of 3 /2 
as in the HM formulation. 
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FIG. 1: Feynman diagrams contributing to the heavy-light (HL) decay constants up to NNLO: (a) and (b) 
are LO and NLO tree level diagrams, loop diagrams (c), (d) and (e) contribute at NLO while diagrams 
(f), (g) and (h) contribute at NNLO. The solid lines denote either HL pseudoscalar or HL vector mesons 
and combinations thereof, dashed lines represent Nambu-Goldstone bosons, the empty (solid) squares and 
empty diamond denote current from the first (third) and second order Lagrangians, and the solid triangles 
denote mass insertions of second chiral order (see Ref. |^]). 

In Eqs. (5,6,7), the superscript in C denotes the chiral order of the corresponding Lagrangian. 
Here we have counted the axial current, the derivative on the NG boson fields, and their masses as 
0(p), as usual. 



To calculate chiral loops, the following LO Lagrangian is introduced tlHd, \2M (in this letter, 
only the relevant terms are explicitly shown): 

£ « = ^Tr[H b H4 abl5 }. (8) 

It describes the interactions between a pair of HL mesons (PP* or p*p*) with a Nambu-Goldstone 
boson <p = k,K,7). In Eq. ®, we have introduced vrtp for the sake of convenience. It should 
be taken as rhn (m B ) at NLO and m D (m B ) at NNLO. In the D meson sector, goD'-w = 9 = 
0.60 ± 0.07 [23], while go*D*-K = 9* is not precisely known. At the chiral order we are working, 
one can take 900*4, = 9dd*-k- If heavy quark spin-flavor symmetry is exact, 933*4, = gB*B*<t> = 
9d*d*4> = 9dd*4>, otherwise deviations are expected. 

The Feynman diagrams contributing to the decay constants up to NNLO 2 are shown in Fig. [TJ 
For the HL pseudoscalar meson decay constants, diagrams (a-g) have been calculated in Ref. [8|]. 
However, diagram (h) that contains two new LECs /Jj and (3 2 was not considered there. Its contri- 



2 The chiral order of a properly renormalized diagram with L loops, Nm (Nh) Nambu-Goldstone boson (HL meson) 
propagators and Vk vertices from fcth-order Lagrangians is n x px — — 2Nm — Nh + J2k ^-Vk- 
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TABLE I: Numerical values of the isospin-averaged masses [24] and decay constants (in units of MeV) 
used in the present study. The eta meson mass is calculated using the Gell-Mann-Okubo mass relation: 
m 2 = {Am 2 K — m 2 )/3. Fq is the average of physical f n , fx and f v . 



m D 



m D 



m B A S (B) A(B) m n m K m v f w 



1972.1 1867.2 102.5 142.6 5331.8 5279.3 88.7 47.5 138.0 495.6 566.7 92.4 1.15/* 



bution to the pseudoscalar decay constant is 



with 



16tt 5 



<j> h {mlml) 



' 2 

mf 



2m t( m l ~ 4m D 



<\> h = 4/3! [m 2 ((4m 2 - m 2 k )B Q {m 2 , m 2 , m 2 ) + A (m 2 )) 
+ (2m 2 -m 2 k )A (m 2 ) 
xB (m 2 , m 2 , m 2 k ) - m • + (im 2 m 2 k + 6mf - 2m 4 k ) 
xA {m 2 ) + 2{5m 2 m 2 k + m$)A (m 2 k ) + m 2 mij , 

where £j can be found in Table 2 of Ref. |8|] with i running over D and D s , j over D* and D*, 
and k over tt, 77, and The functions A = (— 167r 2 )y4 and B = (— 167r 2 )i?o with A and 
B defined in the appendix of Ref. [8]. It should be noted that at NNLO the HL meson masses 
appearing here are the average of the vector and pseudoscalar HL mesons, i.e. m d and m# in Table 
1. For the diagrams contributing to the HL vector meson decay constants, the computation of the 
corresponding diagrams (a, b, e) is the same as in the case of the pseudoscalar decay constants, 
keeping in mind that now a , bo, and 6^ are all understood to be different from those in the 
pseudoscalar sector by heavy-quark spin symmetry breaking corrections. 

The loop diagrams for vector mesons fall into two categories, depending on whether a HL 
vector meson (class I) or a HL pseudoscalar meson (class II) propagates in the loop. For vector 
mesons, the wave function renormalization diagrams (f) yield: 

f i,u _ ( 1 \ f -1 \ d<j)f I ' I \p 2 ,m 2 ,rn 2 k ) 
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3,k 



18F 2 



I67P 



dp 2 
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with 



*\2 



(<f) 



3 (-R 2 + ( m j - m kf){-p 2 i + (rrij + m k f 



xB (p 2 , ml m 2 ) + 3A (m 2 )(-p 2 + m\ - m 2 ) 
-3AK)(P» + ™\ ~ m 2 j) + P 2 (-Pi + 3m 2 + 3m 2 ) 



0/" = 



m 



P „2 



2pf 



3(-2m 2 (p 2 + m 2 ) + (m 2 - p 2 ) 2 + m\) 



xB (p 2 ,m 2 k ,m 2 ) + 3A (m 2 k )(p 2 + m 2 - m 2 ) 
+34 ) (m 2 )(p 2 - m 2 + m 2 ) + 6p 2 (m 2 + m 2 ) - 2p 4 ] , 

where i denotes (D*, D*) and j denotes either (D*, D*) or (£>, D s ). 
Diagrams (g) yields R 9 ' = and 



R 9 ' 



j,k v u 1 



(l67r 2 ) 



<p 9 (m 2 ,m 2 ,m 2 k ) 



with 



^ =-3(( m ,- m ,) 2 - m ?)((m, + m^-m?) 

xB (m 2 , m 2 , m 2 ) + 3A (m 2 )(m 2 + m 2 — ml) 
+3A (m 2 k ){m 2 - m 2 + m 2 .) - 2m? (m? - 3(m 2 + m 2 .)). 



Diagrams (h) give 



with 



if 



a STC ( gmp \ ( ~ l \ Ah'' u ( 2 2n 



m 2 g* 



mp g 



(6m 2 - 3m 2 )A (m 2 ) - 3m 2 m 2 - 2m\ 



+3m 2 k [(4m 2 - m 2 k )B {m 2 , m 2 , m\) + A (m 2 k )] 

40! 3m 2 k {{4m 2 - m 2 k )B {m 2 , m 2 , m 2 .) 
+A (m 2 )) + (6m 2 - 3m 2 )A)(m 2 ) + 6m 2 m 2 



+4m? 



m. 



-6mi(m 2 -4m 2 )B (mlm 2 ,m 2 ) 
+6(3m 2 - m 2 k )(m 2 + m^)A (m 2 ) + 8mfm^ 
+21m 2 mt + 9mf + Q(3m 2 m 2 k + m 4 k )A (m 2 k )^ . 
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As explained in Ref. y8u, mass insertions in diagrams (c, d) generate NNLO contributions. 
Therefore, using mp, a — >■ mp> + A s , mr>* — > rnp> + A, and mo* — > mx> + A + A s for the 
HL meson masses in diagrams (f, g), one obtains the full NNLO results of these diagrams. The 
complete NNLO results for the pseudoscalar and vector HL decay constants are 



ft = a(l + Z l /2) + 5 l + T l + C l + Rl 

F* = ail + (Rf + Rf)/2) + 5 l + tif + R\ + R? + Rf, 

i — I 

where a = a/mp and Z%, Tj, and Cj can be found in Ref. [8]. The "tilde" indicates that one 
has to perform a subtraction to remove the power-counting-breaking terms that are inherent of 



covariant Ch 



see Refs. 



T involving heavy hadrons whose masses do not vanish at the chiral limit (for details 



23]). Furthermore, a second subtraction is needed to ensure that heavy-quark spin- 
flavor symmetry is exact in the limit of infinitely heavy quark masses. Details and consequences 
for phenomenology will be reported in a separate work. After these subtractions the results can 
be expanded in the inverse heavy-light meson mass. In the limit mp — > oo the lowest order 
HMChPT results are recovered. The covariant approach, being fully relativistic, sums all powers 
of contributions in l/m P , which are of higher order in HMChPT. Such a relativistic formulation 
is not only formally appealing. It also converges faster than non-relativistic formulations, such 
as HMChPT and HBChPT. This has been recently demonstrated in the one-baryon sector and in 



2311 . It should be stressed 



heavy-light systems for a number of observables (see, e.g., Refs. [8|, 
that the loop functions are divergent and the infinities have been removed by the standard MS 
procedure, as in Ref. [8|]. 

Now we are in a position to perform numerical studies. We first fix the five LECs, a, bo, b^, 
(5\, and /3 2 , by fitting the HPQCD f D J f D extrapolations |l0|]. The results are shown in Fig. ((2k). 
The NNLO ChPT fits the chiral and continuum extrapolated lattice QCD results remarkably well, 
keeping in mind that the HPQCD extrapolations were obtained using the NLO HMChPT results 
supplemented with higher-order analytical terms [10]. 

In addition to providing the NNLO ChPT results that should be useful for future lattice simu- 
lations of the HL decay constants, a primary aim of the present study is to predict quantitatively 
the SU(3) breaking corrections to f D */fD*, /bJ/b, and /bj//b* from that of the /d s //d- To 
achieve this, one must take into account heavy-quark spin-flavor symmetry breaking corrections 
to the LECs: a, b D , b A , (3 U (3 2 , and gpp,^ (gp+p*^). 

The LEC a is only relevant for the absolute value of the decay constants, therefore it does 



not appear in the SU(3) breaking ratios. However, in the Lagrangian of Eqs. (5,6), one implicitly 
assumes heavy-quark spin symmetry, i.e., d = = 1, which affects the computation of 

loop diagrams (g) for pseudoscalars and (g, h) for vector mesons (see Ref. [8|] for details). Recent 



ll. To be 



quenched LQCD simulations suggest that d is within the range of 1.0 ~ 1.2 II 171 
conservative we allow d to vary within 0.8 ~ 1.2. For b D , bA, Pi, and f3 2 , no LQCD data are 
available. However, the corrections to those constants from heavy-quark spin-flavor symmetry 
breaking are expected to be < 20%. 

The LECs that affect the predicted ratios most prominently turn out to be g and g*, which 
determine the size of chiral loop contributions. In the present case goD*^ is determined by repro- 
ducing the D* meson decay width. Recent rif = 2 LQCD simulations suggest that q R p* T is in 



the range of 0.4 ~ 0.6 Q25 



2711 . We therefore take the central value of 0.516 from Ref. [25] and 



assign a 20% uncertainty. Studies based on QCD sum rules indicate that g and g* could differ by 



10 ~ 20' , |2j,H. We take this into account in our study. 

With heavy-quark spin-flavor symmetry breaking effects on the relevant LECs taken into ac- 
count as described above, we can now make predictions for the ratios of Fd*/Fd*, and 
Fb*/Fb* and their light-quark mass dependencies. The results are shown in Figs. dlbEfcEli). The 
differences between the four ratios are small, at the order of a few percent. Interestingly, the ratios 

1.3 
„ 1-2 
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* HPQCD extrap. 


%, chpT 
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0.2 0.4 0.6 0.8 1 




0.2 0.4 0.6 0.8 1 



0.2 0.4 0.6 0.8 1 



FIG. 2: Light-quark mass evolution of /d s //d> Fd*/Fd*, /bs/Ib, and Fb*/Fb*- The ratio r = m q /m s 
is related to the pseudoscalar meson masses at leading chiral order through = 2Bom s r and m 2 K = 
Bom s (r + 1) with Bq = m 2 /(2m g ), where m s is the physical strange quark mass and m q the average of 
up and down quark masses. The vertical dotted lines denote physical m q /m s . 
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TABLE II: Ratios of /dJ/d, /d*/ fDt±fB s / 7b, and /bj//b* from different approaches. The fD s /fD = 
1.164 from the HPQCD collaboration jlOfl is used as input in our approach. 



Ref. 


fDs 1 Id 






Ib* 1 /b* 


PDG [24J 


1.25(6) 


- 


- 


- 


FCM [30J 


1.24(4) 


1.12 


1.19(3) 


1.15 


RQMQIJ 


1.15 


1.02 


1.15 


1.15 


LFQM [32] 


1.18(1.20) 


1.14(1.18) 


1.24(1.32) 


1.23(1.32) 


QLQCD [17] 


1.10(2) 


1.11(3) 


1.14(3)(1) 


1.17(4)(3) 


QLQCD [18J 


1.11(1)(1) 


1.09(1)(2) 


1.13(1)(1) 


1.14(2)(2) 


LQCD [20J 






1.14(2)(2) 


1.14(2)(2) 


HPQCD [10, 14] 


1.164(11) 




1.226(26) 




NNLO ChPT 


1.17 


1.10(5) 


1.24(4) 


1.20(4) 



of the B meson decay constants are found to be larger than those of their D counterparts, in agree- 
ment with the HPQCD results [Jl3, Jjl. Fully dynamical lattice simulations of the vector meson 



decay constants should provide a stringent test of our predictions. It should be stressed that the 
bands shown in Fig. [2] reflect the estimated effects of heavy-quark spin-flavor symmetry breaking 
from the change of the relevant LECs, in addition to those induced by the covariant formulation 
of ChPT, the use of physical mass splittings and different goD*<j> (9bb*^)- The same is true for the 
uncertainties of our results given in Table II. 



Our predicted ratios at the physical point are compared in Table 



number of other approaches, including the lattice simulations 



113 



18, 



I with the results from a 
2p|l . the relativistic quark 



model (ROM) Bill , the light-front quark model (LFQM) |32|], and the field correlator method 
(FCM) BOIl . 3 Our predictions for the relative magnitude of the fp*/fp* vs. fpj fp ratios agree 
with those of the FCM |30|], the RQM Q3JJ] and LFQM I32n. It should be noted that the results in 
Fig. [2] are obtained with a renormalization scale of 1 GeV I18Q. Uncertainties have been estimated 
changing this scale between \i = mp> and fi = mp for the calculation of D and B decay constants, 
respectively. The changes turn out to be small and are taken into account in the results shown in 



It should be mentioned that the NNLO ChPT predictions cover the NLO predictions within uncertainties. 
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Table EL 

In summary, we have calculated the pseudoscalar and vector decay constants of the B and D 
mesons using a covariant formulation of chiral perturbation theory up to next-to-next-to-leading 
order and found that it can describe well the HPQCD n / = 2 + 1 data on f Ds / f D . Taking into 
account heavy-quark spin-flavor symmetry breaking effects on the relevant LECs, we have made 
predictions for the ratios of /b s //b, /d*//d», and /bj//b* and their light quark mass depen- 
dencies that should be testable in the near future. Our results show that /b s //b > foj fo and 
/dj / Id* < JdJ /d in a large portion of the allowed parameter space. 

This work is supported in part by BMBF, the A. v. Humboldt foundation, the Fundamental 
Research Funds for the Central Universities, the National Natural Science Foundation of China 
(Grant No. 1 1005007), and by the DFG Excellence Cluster "Origin and Structure of the Universe." 
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